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20.  (continued) 

.ritt  gyros  ars  usad  Co  estimate  notion  continuously.  Thraa  Intermittently  comnunlcatlng 
'microcomputers  operate  In  paAaJULtL  to  perforn  the  functions:  (1)  star  image  centroid 
>  determination,  (11)  star  pattern  Identification  and  discrete  attitude  estiaatlon  (sub¬ 
sets  of  nsasured  stars  are  Identified  as  specific  cataloged  stars),  (ill)  optimal 
Kalman  attitude  motion  estimation/integration.  The  system  is  designed  to  be  self- 
calibrating  with  provision  for  routine  updating  of  Interlock  angles,  gyro  bias  parameters, 
and  other  system  calibration  parameters.  For  redundancy  and  improved  precision,  two 
optical  ports  are  employed.  This  interim  report  documents  Phase  I  of  a  three  phase 
effort  to  research,  develop,  and  laboratory  test  the  basic  concepts  of  this  new  system. 
Included  la  Phase  1  is  definition,  formulation,  and  test  of  the  basic  algorithms, 
including  preliminary  implementations  and  results  from  a  laboratory  microcomputer 
system. 
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SUMMARY 


The  primary  result*  of  the  Phase  1  effort  are  the  following: 

1.  Development  of  an  approach  for  real  time  on-board  estimation 
of  spacecraft  orientation  with  aub  five  arc-second  precision. 

2.  Detailed  formulation,  testa  and  validations  of  an  efficient 
and  reliable  star  pattern  recognition  strategy  appropriate  for 
use  with  charged-coupled-device  (CCD)  array- type  a tar  sensors. 

3.  Formulation  and  preliminary  tests  of  a  motion  integration/ 
Kalman  filter  algorithm  to  integrate  gyro  measured  angular 
rates  and  (by  sequential  processing  of  the  discrete  orien¬ 
tation  information  available  from  the  star  sensing,  identi¬ 
fication,  and  attitude  determination  process)  provide  optimal 
real  time  estimates  of  spacecraft  orientation  and  angular 
velocity. 

4.  Definition  of  software  and  hardware  for  a  laboratory  micro¬ 
computer-system-based  study  to  permit  realistic  ground-based 
validation  and  optimisation  of  the  algorithms  and  computer 
configuration. 

These  results  are  discussed  in  detail  herein. 


Ill 


1.0  INTRODUCTION  AND  SYSTEM  OVERVIEW 

The  purpose  of  Che  research  activities  reported  upon  herein  is  to 
research  and  develop  a  system  which  has  potentially  profound  impact  upon 
satellite-based  technologies  requiring  precision  knowledge  (sub-five- 
arc-  second  uncertainty)  of  spacecraft  attitude. 

The  focus  of  the  research  has  been  the  optimal  exploitation  of 

array-type  Charged  Coupled  Devices  (CCD)  as  the  star  light  sensitive 

element  in  a  new  generation  of  star  sensing  and  star  pattern  recognition 

systems;  these  systems  hold  the  promise  of  autonomous,  on-board,  near- 

real-time  determination  of  spacecraft  attitude,  with  state-of-the-art 
1-3 

precision  .  The  heart  of  the  measurement  hardware  designs  (several 
similar  configurations  are  currently  under  development  and  test)  is  a 
CCD  array  or  matrix  of  about  200,000  light-sensitive  elements,  accurately 
imbedded  (to  1  part  in  10,000)  in  a  microcircuit  chip.  The  CCD  array 
serves  as  the  "film"  of  a  new  generation  of  "digital  stellar  cameras" 
which  can  accurately  digitize  6tar  Images  with  negligible  lag  from  real 
time.  The  digitized  starlight  is  analyzed  by  a  programmable  microprocessor 
which  extracts  star  image  coordinates  and  magnitudes  and  associates  a 
time  with  these  measurements.  The  microprocessor  outputs  star  image  and 
magnitude  data  in  near-real  time  (current  projections  indicate  5  or  more 
star  "photographs"  per  minute).  The  measured  star  patterns  output  from 
the  microprocessor  can  be  telemetered  to  ground,  but  ultimately,  it 
appears  feasible  to  use  methods  described  herein  to  determine  spacecraft 
orientation  from  these  data  with  sub-five-arc-second  precision,  within 
the  constraints  of  on-board  computation,  and  thereby  open  the  door  for  a 
truly  dramatic  breakthrough. 
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With  reference  to  Figures  1.1,  1.2,  and  1.3,  ve  discuss  a  modular 
parallel  process  approach  to  the  attitude  determination  problem.  There 
are  four  major  sub-systems,  viz: 

(1)  Two  CCD  matrix  arrays  and  associated  electro-optics  for  sensing 
and  digitizing  starlight.  Two  identical  optical  ports  having 
effective  focal  lengths  of  about  70  mm  focus  starlight  from  7° 
x  9°  fields  of  view  onto  the  CCD  arrays.  The  fields  of  view 
can  be  separated  by  any  convenient  angle,  but  we  have  nominally 
used  90°.  The  CCD  arrays  are  taken  as  Fairchild's  11.4  mm  x 
8.8  mm  matrix  arrays  containing  a  grid  of  488  x  380  silcon 
pixels.  The  starlight  is  slightly  defocused  on  the  CCD  image 
planes  (to  facilitate  more  accurate  "centroiding"  of  image 
locations);  typical  stars  illuminate  16  to  25  pixels.  The 
image  processing  takes  place  in  two  stages.  The  first  stage 

is  a  high  speed  analog  scan  to  identify  the  pixel  addresses 
and  A/D  convert  the  response  levels  which  are  greater  than  a 
preset  analog  level.  The  analog  scans  for  the  two  fields  of 
view  are  controlled  by  a  common  clock  to  make  the  digitized 
data  valid  at  the  same  instant. 

(2)  Microprocessor  A  performs  the  program  Process  A,  which  is  the 
digital  portion  of  the  image  process  function.  Process  A  is  a 
sequence  of  calculations /logical  decisions  requiring  as  input 
the  A/D  converted  response  levels,  pixed  addresses,  and  time 
from  the  CCD  response.  The  output  data  of  Process  A  are  the 
Interpolated  image  centroids  and  Instrument  (relative)  magni¬ 
tudes  of  each  valid  star  image.  An  average  of  about  six  stars 
will  be  detected  on  each  CCD  array.  Upon  ouput  from  Process 
A,  the  star  image  centroids  will  have  been  enhanced  by  all 
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preflight  calibration  data,  and  any  easily  detectable  spurious 
images  will  have  been  deleted.  The  time  recorded  when  the 
data  were  scanned  from  the  CCD  arrays  remains  associated  with 
data  and  is  transmitted  along  with  the  image  centroids  and 
magnitudes  to  Process  A's  output  queue.  The  output  queue  of 
Process  A  will  be  continuously  accessible  by  Process  B  and/or 
telemetry  systems;  the  information  will  be  updated  (over¬ 
written)  at  a  to-be-determined  interval  (about  five  times  per 
minute).  The  data  compression  affected  by  immediate  extraction 
of  image  centroids  and  magnitudes  from  the  pixel  response 
(and,  by  virtue  of  the  high  speed  analog  scan,  considering 
only  the  significantly  illuminated  pixels)  is  most  important 
to  the  concept  of  an  autonomous  on  board  system;  the  trans¬ 
mission:,  storage,  and  subsequent  analysis  of  the  large  volume 
of  pixel  data  is  thereby  made  unnecessary. 

(3)  Microprocessor  B  is  programmed  to  execute  Process  B  which 
accepts  as  inputs: 

•apriori  attitude  estimates, 

*apriori  estimation  of  interlock  angles  between  the  two 
image  planes 

■star  image  coordinates  and  magnitude  data,  usually 
one  data  set  from  each  CCD  array,  from  Process  A  (port 

ID. 

Through  a  sequence  of  calculations/decisions  (including 
systematic  access  of  a  star  catalog  containing  the  direction 
cosines  and  magnitudes  of  the  5000  brightest  stars).  Process  B 
performs  two  primary  functions: 

‘Identifies  the  measured  stars  in  each  field  of  view-of- 
vlew  as  specific  cataloged  stars  (using  a  "star  pair 
correlation"  technique),  and 
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•Determines  the  spacecraft  orientation  angles  and  field- 
of-viev  interlock  angles  (which  causes  the  simulated 
images  of  identified  stars  to  overlay  the  measured 
images  in  the  least  square  sense). 

The  spacecraft  orientation  output  from  Process  B  will  be 
available  continuously,  and  will  nominally  be  updated  two  or 
more  times  per  minute.  As  soon  as  the  attitude  determination 
process  is  completed,  the  previous  attitude  angles  and  associ¬ 
ated  time  are  overwritten  by  the  new  attitude  angles  and  times 
(at  data  port  III,  at  a  nominal  rate  of  2  or  more  readouts  per 
minute) . 

(4)  Microprocessor  C  is  programmed  to  accept  the  discretely  valid 
attitude  determinations  of  Process  B  and  A/D  converted  gyro 
measurements  of  angular  velocity  to  determine  a  best  real  time 
attitude  and  angular  velocity  estimates  for  possible  feedback 
to  the  control  system,  for  feedback  to  Process  B,  and  for 
telemetry  to  ground. 

In  Phase  I  (just  completed)  of  this  research  and  development  effort, 
have  concentrated  upon  the  following  tasks: 

(1)  Development  and  test  of  methods  and  algorithms  for  Process  B, 
for  star  pattern  recognition  and  attitude  determination. 

(2)  Development  and  test  of  methods  and  algorithms  for  Process  C, 
for  continuous  real  time  attitude  estimation  from  Process  B 
discrete  attitude  estimates  and  rate  gyro  measurements. 

(3)  Evaluating  the  designs  and  error  characteristics  of  competing 
CCD  star  tracker  designs,  in  order  to  develop  realistic  CCD 
output  simulations  for  input  to  Procsss  A  and  B. 
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(4)  Simulation  of  instrument  response  characteristics  on  the  CCD 
array,  centroiding  errors,  etc. 

(5)  Creation  of  mission  star  catalogs. 

(6)  Identification  of  requirements  for  a  microcomputer  system  to 
conduct  tests  of  Processes  B  and  C  within  the  precision,  speed, 
and  storage  constraints  of  an  on-board  computer. 

To  date,  the  overwhelming  majority  of  the  computational  studies  have  been 
carried  out  on  the  VPI  &  SU  IBM  370/158  computer  system.  However,  we 
have  recently  begun  conversions  and  report  some  limited  results  of  studies 
done  on  a  HP  9845S  microcomputer  system.  During  Phase  II,  both  Processes 
B  and  C  will  be  implemented  on  this  computer. 

In  the  text  of  this  report,  we  discuss  the  above  tasks  and  their  results, 
status,  and  outlook  for  Phases  II  and  III. 
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UVASTAR  An  electro-optical/ software  system  capable 
of  raal  tine  readout  of  digitized  atar  coordinates, 
and  ultisatal7,  autonomous,  near-real  time  star 
pattern  recognition  and  attitude  determination. 
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FIGURE  1.2:  Format  Ion  of  I  mat*  on  Cho  CCD  Array 


2.0  COORDINATE  FRAME  NOMENCLATURE  AND  KINEMAT1CAL  RELATIONSHIPS 

Important  to  any  complex  angular  motion  problem  le  the  choice  of 
orientation  parametera  and  a  clear,  ayatematlc  nomenclature  convention. 

The  choice  of  orientation  parametera  la  critical  in  that  numerical 
difficulties  may  be  avoided  entirely  or  at  least  relegated  to  attitudes 
which  will  never  occur  in  the  problem  at  hand;  a  self-evident  nomenclature 
is  useful  (particularly  when  several  coordinate  systems  are  defined)  not 
only  to  readers  but  to  the  investigators  themselves. 

Throughout  Phase  I  of  this  research  effort,  Euler  angles  have  been 
utilized  for  coordinate  frame  orientation.  It  is  a  well-known  fact  that 
each  of  the  twelve  possible  Euler  angle  sets  possesses  a  kinematlcal 
singularity.  In  order  to  render  the  singularity  problem  as  harmless  as 
possible,  the  Euler  angle  sequence  for  a  given  coordinate  frame  relative 
orientation  has  been  selected  so  that  encountering  the  singular  orientation 
is  highly  unlikely.  Future  efforts  may  involve  the  uBe  of  a  four  parameter 
orientation  description  (such  as  Euler  parameters)  so  as  to  avoid  the 
singularity  problem  entirely. 

Unless  there  are  reasons  to  the  contrary,  the  1-2-3  Euler  angle 
sequence  will  be  preferred.  For  future  reference,  the  3x3  transformation 
matrix  [DE]  orienting  a  frame  D  characterized  by  the  orthonormal  vector  set 
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is  compactly  written  as  the  matrix  product 

U>  -  tOXW^Hfs)- 

Obtainable  from  an  elementary  kinematic  analysis  of  the  frames  displayed 
in  Fig.  2.1  is  the  following  matrix  expression  for  the  components  (along 
the  (d)  unit  vectors)  of  the  angular  velocity  of  frame  D  with  respect  to 
frame  E: 

D/E. 


D/E 


80 


de 


:C*d« 

8^de 

o' 

»*d. 

c*de 

0 

0 

1 

(2.5) 


Finally,  both  members  of  Eq.  (2.5)  may  be  premultiplied  by  the  inverse  of 
the  3x3  matrix  on  the  right  hand  side  to  obtain  expressions  for  the  time 
rates  of  change  of  the  Euler  angles : 
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D/E1 


Contained  in  the  above  discussion  are  sub-and  superscript  conventions 
which,  while  appearing  unnecessary  in  this  simple  context,  will  be  utilized  to 
great  advantage  in  conjunction  with  the  several  coordinate  frames  about  to 
be  introduced. 

We  now  turn  to  the  first  of  the  specific  coordinate  frames  selected  for 
use  in  the  present  problem.  The  inertial  frame  N  (see  Fig.  2.2)  la  non-rotating 
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and  geocentric;  ic  forms  the  beats  for  the  on-board  star  catalog  and  is 
characterised  by  the  unit  vector  set 
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An  "orbit  inertial"  frame  O'  (see  Fig.  2.2)  is  defined  by  the  unit  vectors 
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where  the  subscript  o  denotes  Initial  conditions.  The  {£* }  unit  vectors 
are  projected  onto  the  {n}  unit  vectors  by  the  direction  cosine  matrix 
[O'N]  as 

(o')  -  [O'N] (n)  (2.9) 

The  nine  constant  elements  of  [O'N]  are  calculated  from  the  orbit  state 
vector  by 
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The  orbit  frame  0  (see  Fig.  2.3)  rotates  about  £2»  the  orbit  normal, 
with  the  constant  angular  velocity  (relative  to  N  or  O') 

(2.12) 
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CM  is  the  earth  gravitational-mass  constant  and  a  is  the  semi-major  axis 

of  the  orbit.  The  (o)  unit  vectors  are  written  in  terms  of  the  (o')  unit 

* 

vectors  by 
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(2.14)  yields 
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(2.16) 


The  next  coordinate  frame  (G)  considered  is  one  attached  to  the  strapped- 
down  gyroscopic  assembly  and  hence  body-fixed  but  not  necessarily  principal. 
The  G  frame  is  generally  considered  to  be  the  primary  spacecraft  axes.  The 
unit  vector  set  <£}  of  Fig.  2.3  is  characterised  by  £A  being  the  roll  axis 
(nominally  along  the  velocity  vector),  £2  being  the  pitch  axis  (nominally 
along  the  orbit  normal),  and  ^  being  the  yaw  axis  (nominally  along  the 

radius  vector).  The  1-2-3  angles  4  ,  8 _ ,  and  *  orienting  G  relative  to 

go  go  go 

0  are  banco  normally  small,  thereby  avoiding  the  kinanatical  singularity  at 
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biiww.- 


6  ■  90*.  The  {g}  unit  vectors  ere  projected  onto  the  {o}  unit  vectors 

go 

by 

(*>  -  i<»<vVV1{£l*  (2,17) 

Utilising  Eq.  (2.16): 

{*}  -  [GO] [00' ] [O'N] {n} .  (2.18) 

Ue  now  introduce  the  frames  A  and  B  associated  with  the  star  sensor 

fields-of-view  A  and  B.  Frame  A  (see  Fig.  2.4)  is  characterized  by  a^  being 

along  the  fleld-of-view  A  boresight,  with  a^  and  a^  being  image  plane 

axes.  The  angles  $  ,  8  ,  and  orient  the  {a}  unit  vectors  relative 

ag  ag  ag  — 


to  {£}  via  the  3-1-3  rotation 

{a}  -  lAG(*ag,eag.*#g)HA>  (2.19) 

-  [AG] [GO] (o)  (2.20) 

-  [AG]  [GO]  [00' Ho'}  (2.21) 

-  [AG] [GO] [00' ] [O'N] {n} .  (2.22) 


A  3-1-3  sequence  has  been  chosen  for  thi6  orientation  interlock,  with  the 
nominal  values 


(♦ag»eag’*ag)  "  <90°,135#,0») .  (2.23) 

Frame  B,  for  purpose  of  recovering  the  interlock  angle  between 
fields-of-view  A  and  B,  is  oriented  relative  to  frame  A  by  the  3-1-3  angles 

**d  *b. vl* 

(b)  -  [BA(*ba,eba,*ba)]{a}.  <2.2A) 

For  our  study,  the  interlock  angles  have  been  assigned  the  nominal  valuer 
(see  Fig.  2.4). 

(*ba’eba»*ba)  "  ("90°’  90°’  90#)*  (2.25) 

The  orientations  prescribed  by  Eqs.  (2.23)  and  (2.25)  correspond  to  a^  (the 

field-of-view  A  boresight)  being  "yawed"  about  ^  by  45*  from  and  b^ 
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(the  f ield-of-view  B  boresight)  being  "yawed"  about  by  -45°. 

Finally,  we  introduce  an  "orbit  A  frame",  designated  A’,  which  is 

coincident  with  frame  A  at  the  initial  time  t  ,  but  which  thereafter 

o 

rotates  about  the  orbit  normal  "  o^  with  the  constant  angular  velocity 

A'/N  A' /O'  o  «  n  ,, 

w  ■  u  “  -  iio^.  (2. 

Frame  A*  represents  the  nominal  position  of  frame  A  and  is  oriented 

relative  to  frame  0  via 

{a’}  -  [A'0(*  .  ,0  .  .  )Ho},  (2.27) 

—  a  o  a  o  a  o  — 

Where  the  1-2-3  Euler  angles  ♦a,o,  ea'0»  anc^  ^a'0  are  constants.  The 
relationship  between  the  frames  A  and  A'  is  given  by 

{a}  -  (2. 


(2.26) 


(2.28) 


wherein  the  three  1-2-3  angles  $  ,,  0  ,,  and  ^  ,  are  normally  small. 

33  33  33 

Table  1  summarizes  the  coordinate  system  definitions  and  relationships 
discussed  thus  far. 

The  interlock  angles  ($  ,6  ,  ip  )  and  (4>,  ,  6.  ,  tk  )  will  vary 

ag  ag  ag  ba  ba  ba 

about  their  nominal  calibrated  values  due  to  thermal  cycling  or  any  other 

structural  deformation.  The  maximum  amplitude  of  such  variations  is 

believed  to  be  a  few  arc-minutes.  As  is  evident  in  the  developments  of 

Sections  3  and  4,  it  is  possible  to  update  the  estimates  of  the  interlock 

angles  ($.  ,  0.  ,  <|».  )  between  the  star  sensor  coordinate  systems, 

ba  ba  ba 

However,  we  have  not  yet  structured  the  estimation  algorithms  of  Process 

C  to  provide  correction  of  the  interlock  angles  ($  ,0  ,  4’  ).  It  is 

ag  ag  ag 

believed  that  noise  in  typical  rate  data  (around  one  arc-sec/sec)  is  too 
large  to  allow  significant  refinement  of  ($  ,6  ,  )  unless  Processes 

B  and  C  can  be  implemented  in  such  a  fashion  that  several  Process  B 
outputs  per  minute  can  be  passed  through  the  Kalman  filter  of  Process  C. 

As  is  noted  in  Section  4,  these  interlock  angles  are  correlated  with  the 


estimated  rate  bias  parameters,  a  fact  which  also  hampers  accurate 


estimation 
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Useful  for  consideration  in  Process  C  (see  Section  4)  are  the  time 


rates  of  change  of  various  Euler  angle  sets  in  terms  of  the  Inertial 
angular  velocity  components  of  frame  G.  Specifically,  we  will  develop 
expressions  for 


(1)  (♦  ,6  ,*  ) 

gn  gn  gn 

(2)  (L,0  ,i  ) 

an  an  an 

(3)  (*  .0  ) 

go’  go,rgo' 

(4)  (♦aa..®aa«.*aai>* 

The  relationships  for  the  time  rates  of  change  of  the  angles 

6  ,  6  ,  and  <|/  follow  directly  from  the  form  of  Eqs.  (2.6): 

vgn*  gn’  ygn  7 


i  $ 

l  gn 

1  e 

(0 

1  gn 

1 

C0 

I  e 

gn 

\i(< 

T  nn' 

gn 


'c<(i 


gn 


-8l(l 


gn 


C0  811/  C0  Clil 
gn  gn  gn  ygn 


0 

0 


80  Cl)/ 
gn  gn 


gn  gn 


gn 


■e 

W1 

G/N 

“2 

G/N 

U>-j 

(B(0  ,1)/  )]  (w 

gn  gn  ' 


G/N 

1 

G/N 

2 

G/N 


Note  that  0  ,  for  the  case  of  a  near  circular  orbit  and  one  vehicle 

gn’ 

rotation  about  Per  orbit,  will  go  to  90°  twice  per  orbit,  causing 

CO80  to  vanish  with  a  resulting  singularity  in  Eqs.  (2.29). 
gn 

The  relationships  for  the  time  rates  of  change  of  the  angles  <p 

(2 

G/N 


6  ,  and  t|/  can  also  be  obtained  from  Eq.  (2.6) 

an  an 


an 


an  i 


anl 


-  [  B  (  0  ,i(i  ) )  [  AG  (  $  ,0  ,i|<  ))  <  to  „ 

1  v  an*yan  11  ag’  ag,rag  \  2 


an 


iw- 


G/N ! 
G/J 


(2. 


(2.30) 


The  factor  of  (AG)  on  the  right  side  of  Eq.  (2.30)  is  necessary  to  convert 
the  from  components  along  (jj)  axes  to  components  along  (a)  axes.  Due 

to  the  form  of  the  matrix  B  (defined  in  Eq.  (2.6)  and  repeated  in  Eq. 
(2.29),  singularities  would  occur  twice  per  orbit  for  a  vehicle  which 
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29) 


rotates  about  once  per  orbit.  Hence,  although  fairly  simple  in  form, 

Eqs.  (2.29)  and  (2.30)  share  similar  difficulties. 

We  now  proceed  to  develop  equations  analogous  to  (2.29)  and  (2.30) 

for  the  angles  (t  ,0  ,  i 1>  )  (Case  3) .  The  angular  velocity  of  G 

go  go  go 

relative  to  N  may  be  written  as 

uG/N  -  wG/O  +  uiO/N  (2.31) 


or 


hi,  G/N  +  ui_G/N  +  ui.G/N 
1  ii  2  &2  3 


o,(1)+0  0_(2)  +  t  +  «o,,  (2.32) 

gcr-1  go-2  r  go  -2* 


H  *2  '  *3 

where  the  parenthetical  superscripts  indicate  intermediate  frames  in  the 

0  to  G  transformation.  Expressing  o  ^  ^ ,  Oj^,  and  o^  in  terms  of 


,  and  enables  one  to  equate  coefficients  of  like  unit  vectors, 
resulting  in  three  scalar  equations  which  have  the  matrix  form 


c*  sili  +  st  s0  ciji 
1  Tgo  Tgo  Tgo  go  go 


11  \  c<t>  Clfl  -  St  S0  St  , 
'  go  go  *go  go  rgo| 

-St  C0 
go  go 


(2.33) 


For  the  case  of  one  vehicle  rotation  about  Per  near  circular  orbit, 

the  angles  t  ,  6  ,  and  t  remain  small  so  that  Eqs.  (2.33)  are  non- 
6  Tgo’  go’  go 

singular  and  nearly  linear  (important  in  the  implementation  of  estimation 
algorithms) . 

•  a  • 

To  develop  the  Case  4  equations  for  t  ,,  0  , ,  and  t  , ,  we  begin 

8A  AA  Afl 

with 


G/N  A/A'  .  A'/N  . 

ui  -  hi  +  w  .  (2.34) 

Proceeding  in  a  manner  analogous  to  that  outlined  for  Case  3,  we  arrive  at 


*,..11 


[AG(t  ,0  ,t  ) 
ag’  ag  ag 


21 


FOV  B 


»*H  .  L„  ♦  BMU  .  !■„  ♦  B»13  ■  L,, 
b  B,,31  '  L1J  +  BN32  '  L2j  +  B33  •  ■•3]  “pb 


(3.2a) 


'  ,U"Ctl°n  \a WV  V-V 

orientation  interlock  4.' s  dir.  cosines  FOV  B 
^  's  of  FOV  A  from  B  to  A  toward  calibration 

jth  star  constants 


y  .  f  BM21  ■  L1J  +  BM22  '  L21 
bJ  b  >»»31  .  L;j  +  BN32  .  L 


l21  +  BN23  •  L3i 

L2j  +  BN33  .  L3j  J  +  ypb 


(3.2b) 


-  function  <*a^an**an;  ♦bAA5  L1J*V  L3j :  ypb’fb> 

orientation  interlock  dir.  cosines  FOV  B 

^  's  of  ^ 's  from  toward  jth  calibration 

FOV  A  B  to  A  star  constants 

[BN]  -  *  '“I1**" 


|BA<»b.-Bb.-*ba>l 


'  c*ba  •♦b« 
-♦b.  c*b. 


C0. 

ba 

80. 

ba 

-sQ. 

ba 

c 0.  . 
ba 

ba 

8*ba 

0" 

ba 

C^ba 

0 

(3. 2c) 

! 

0 

1. 

in 

FOV  B. 

In  trying  to  compare  measured  image  centroid  coordinates  from  Process  A 
with  computed  image  coordinates  from  Eqs.  (3.1)  and  (3.2),  one  is 
immediately  faced  with  several  difficulties: 


74 


’  C(^  S«P  +  8$  80  Clli 

80  80  *80  80  *80 I 


-  ntA'0(^,o,ea,o^a,o)]  {c^aoc*ao  -  s^se^ 


801 


(2. 


-84  c0 
80  80 


Like  Eqs.  (2.33),  Eqs.  (2.35)  are  more  complex  than  Eqs.  (2.29)  or  (2.30) 
but  offer  the  advantages  of  near  linearity  and  no  singularities  for  the 
case  of  one  vehicle  rotation  about  Per  orbit. 


35) 
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3.0  STAR  PATTERN  RECOGNITION/ATTITUDE  DETERMINATION 
3.1  Some  Preliminary  Remarks 

Figure  3.1  displays  a  typical  set  of  measured  image  coordinates, 
a  simulated  set  of  coordinates  (based  upon  apriori  angle  estimates)  and 
an  uncertainty  region  on  the  celestial  sphere. 

The  cataloged  star  images  are  projected  from  the  celestial  sphere 
onto  the  CCD  image  planes  via  the  colinearity  equations 


FOV  A 


.  ANU  .  Lu  +AN12  .  Ln  +  ANn  ■  L 
Xai  alAN31  .  +  AN32  .  L21  +  AN33  .  L3iJ  Xpa 


(3.1a) 


"  funCtion  Lli’L2i,L3i’ 

orientation  dir.  cosines  FOV  A 
■dL’s  of  FOVA  to  ith  star  calibration 

constants 


f  AN21  ‘  Lli  *  AN22  ‘  L2i  *  AN23  *  L3i. 
a‘AN31  .  Lu  +  AN32  .  L2i  +  AN33  .  L3iJ  ypa 


(3.1b) 


-  function  (*an.ean,*an;  Lu.  L^.  L^J  ypa.fa) 


orientation  star  dir. 
of  cosines 

FOV  A 


FOV  A 

calibration 

constants 


[AN]  -  l AM (♦  . 0  . ♦  ) 1  - 

an  an  an 


ciji  a\p  01  fc6  0  -s0 
an  an  II  an  an 


■si p  cip  0(|0  10  I  1 0  c<p  sip 

an  an  II  II  an  an 

0  0  lJ  LsO  0  cO  J  LO  — s4>  c<J> 

an  an  an  an 


0  0 


(3.1c) 


cosa^  cosfi^ 
sina^  cos6^ 


*  direction  cosines  toward  the  ith  star 
appearing  ln  FOV  A,  referenced  to  the 
Inertial  axes  of  the  mission  star  catalog 
(see  appendix  1) 


(3. Id) 
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(1)  The  orientation  angles  (♦  ,6_  ,*  ),  In  general,  are  not 

an  an  an 

known  with  enough  precision  for  the  projected  catalog  stars  to 
overlay  measured  images. 

(2)  The  catalogued  stars  corresponding  to  measured  images  will 
not  generally  be  known  beforehand,  thus,  the  particular  star 
direction  cosines  which  should  be  used  in  (3.1)  and  (3.2)  to 
generate  coordinates  for  comparison  with  particular  measured 
values  will  not  be  known  apriori. 

One  brute-force  approach  would  involve  many  trial  and  error  iterations 

of  ($  ,6  , )  and,  for  each  orientation  trial,  one  could  map  all 

an  an  an 

stars  in  a  particular  region  of  the  celestial  sphere  through  Eqs.  (3.1) 
and  (3.2)  and  then  attempt  to  correlate  them  with  the  measured  pattern. 
Such  an  approach  would  be  ill-advised  and  completely  unsatisfactory  in 
view  of  its  inefficiency,  especially  for  on-bcard  computation. 

The  key  to  efficient  star  identification  is  to  take  advantage  of 
the  sub-ten-arc  second  precision  of  Process  A;  the  angles  between  pairs 
of  measured  stars  are  very  well  determined  by  the  measured  coordinates 
and  can  be  used  to  identify  the  corresponding  catalogued  stars.  The 
cosine  of  the  angle  between  a  typical  pair  of  measured  stars  can  be 
computed  from  the  measured  image  coordinates  as 


C^  =  COS0 


(x±x  +  y±y  +  f  ) 


1;i  /  2~  2  "  ,2W  2~  2  ,  J 

/<xi  +  yi  +  f  )<xj  +  +  f 


(3.3) 


The  cosine  of  the  angle  between  a  typical  pair  of  cataloged  stars  can  be 
computed  from  the  cataloged  data  as 

CIJ  ■  •  Vl,  +  L2,L2,  +  Vv  (3-4) 

X  J  i.  J  1  J 

The  pattern  recognition  logic  developed  below  makes  use  of  the  smallness  of 


the  difference  between  (3.3)  end  (3.4)  es  e  Men*  to  tentatively  identify 
the  measured  etare  in  the  catalog. 

The  sequence  of  calculations  and  logical  decisions  for  accomplishing 
the  pattern  matching  and  attitude  angle  estimation  is  not  unique,  but  we  feel 
an  autonomous  on-board  algorithm  should  have  the  following  features: 

(1)  Allow  a  wide  variation  in  the  quality  of  aprlori  orientation 
uncertainly  (from  say,  several  degrees  (start  up  situation) 

to  several  arc  seconds  (steady  state  mode)],  yet,  he  relatively 
efficient  and  highly  reliable  in  each  case. 

(2)  Allow  a  "reasonable"  mis-match  between  the  set  of  cataloged 
and  measured  stars  (  to  reflect  the  realistic  occurence  of 
occasional  spurious  measured  images,  measurement  of  binary 
stars  (which  we  have  deleted  from  the  mission  catalog  to 
avoid  the  problems  associated  with  them),  and  uncertainties 

in  the  spectral  properties  of  the  CCD  array  and  optical  system 
and  the  cataloged  information]. 

(3)  Permit  continuous  re-calibration  of  the  interlock  angles 
between  FOV  A  and  FOV  B,  to  compensate  for  the  common  event  of 
significant  thermal  deformation  cycling. 

(4)  Seemingly  in  conflict  with  the  first  three  requirements,  the 
steady-state  "cycle  time"  (to  complete  the  transformation  from 
image  coordinate  measurements  to  output  orientation  and  inter¬ 
lock  angle  measurements)  should  he  kept  small;  it  is  felt 
most  desirable  that  the  cycle  time  be  less  than  30  seconds 
(within  the  constraints  of  on-board  computation). 

A  macroscopic  "functional  logic  flow"  of  the  present  Process  B 
software  (presently  operational  on  the  IBM  370/158,  computer  to  be 
converted  to  the  HP  9845  microcomputer  during  Phase  11)  is  shown  in  Figure 
3.2,  This  system  has  been  Implemented  in  FORTRAN  and  (on  the  basis  of 
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simulations,  as  suunarised  in  Saction  5.0)  appears  to  meet  the  first  three 
of  the  above  objectives  quite  satisfactorily.  The  degree  to  which 
objective  (4)  can  be  net  will  be  studied  in  detail,  by  implementing 
and  optimising  this  software  on  the  HP  9845S  computer  during  Phase  11. 
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FIGURE  3.2  PROCESS  B 

STAR  IDENTIFICATION/DISCRETE  ATTITUDE  ESTIMATION 
FUNCTIONAL  LOGIC  FLOW 


Begin  New  Cycle 

©i 

Input  FOV(A)  Image  Centroid  Measurements  from  Process  A 
Input  Apriori  Attitude  Estimates  from  Process  C  of  FOV(A) 

-4=  s  ($,e,40. 


®1 


an 


Retrieve  Cataloged  Stars  within  an  uncertainty  border  of 
FOV(A)  from  mission  catalog. 

@-j 

Sort  Stars  by  their  angle  off  the  estimated  boresight  (small  to  large). 

•1 

Compare  Sequence  of  Cosines  of  the  angles  between  pairs  of  cataloged 
stars  to  the  cosines  of  angles  between  pairs  of  measured  stars. 

A  failure  to  find  a  match  returns  control  to  step  Q  . 


Least-Squares  Differential  Correction  (LSDC)  to  determine  ($,6,$) 
which  yields  the  best  overlay  of  simulated  image  coordinates  anda 
measured  coordinates  for  the  identified  stars.  Check  for  additional 
overlays  of  cataloged  and  measured  stars  to  confirm  orientation.  If 
not  confirmed,  continue  with  (D  • 


©-J 


Retrieve  Cataloged  Stars  within  uncertainty  border  of  FOV(B)  from 
mission  catalog.  Use  ($,9, 41)  from  step  and  apriori  interlock 

angles  ($,0,$).  toaorient  FOV(B)  for  this  step. 


Repeat  steps  Q  _-  ©  for  FOV(B)  holding  ($,0,<|O.  constant. 


LSDC  to  determine  ($,9,$)  which  yields  best  overlay  for  stars  in 
both  FOV(A)  and  FOV(B)  simultaneously.  Use  up  to  5  confirmed  stars  from 
each  FOV  and  hold  ($,9,$)^  constant. 


®1 


LSDC  to  determine  ($,9,<0.  which  yields  best  overlay  for  stars  in 

FOV(B)  only.*  Hold  ($,0,$)  constant. 

an 


©-j 


Output  (to  Process  C  buffer)  the  best  estimates  of  ($,9»40  ,  ($,0,40. 

and  corresponding  covariance  matrices.  *n  * 


®1 

Return  to  step  for  new  data. 
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4.0  ATTITUDE  PREDICTION/OPTIMAL  SEQUENTIAL  ATTITUDE  ESTIMATION 

4.1  Preliminary  Comments  on  Definition  of  the  State  Vector 

Process  B,  as  described  in  section  3.0,  provides  discrete  estimates 

of  orientation  angles  <4  ,  0  ,  )  orienting  FOV  A  in  inertial  space 

ax\  An  an 

and  interlock  angles  G^a ,  ij/^)  between  FOV  B  and  FOV  A.  These 

estimates  are  valid  at  instants  which  lag  real  time  by  an  amount  equal  to 

the  real  elapsed  time  required  to  perform  the  computations  of  Process  B. 

During  this  finite  time  interval  (typically,  a  fraction  of  a  minute), 

the  gyro  measurements  of  angular  rates  can  be  integrated  to  estimate 

vehicle  orientation.  Unfortunately  several  sources  of  uncertainty  corrupt 

these  data  in  real  life  applications: 

rate  biases,  noise,  A/D  conversion  errors,  time  delays,  and  lack 
of  orthogonality  (gyro  instrument  errors) 

uncertain  interlock  angles  between  the  gyro  axes  and  the  FOV  A 
CCD  image  plane  axes. 

Of  these  error  sources,  the  most  significant  are  the  rate  biases. 
Fortunately,  the  gyro  rate  measurements  will  typically  be  Integrated  for 
only  a  fraction  of  a  minute  before  receipt  of  heavily  weighted  angles  from 
Process  B,  so  small  uncompensated  gyro  errors  will  not  be  reflected  in  a 
long-term  degradation  of  precision.  Any  attitude  estimation  algorithm 
which  ignores  one  or  more  of  the  above  problems  (for  the  purpose  of 
establishing  a  realizable  algorithm)  should,  however,  be  evaluated  by 
establishing  "truth  models"  in  which  realistic  rate  data  can  be  processed 
by  a  sub-optimal  Process  C  algorithm  then  comparing  the  estimated 
orientation  history  with  the  truth. 

If  on-board  computational  constraints  were  not  present,  one  would 
probably  consider  developing  an  algorithm  for  estimating  the  relatively 
high  dimensioned  state  vector 
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orientation  biases  in  interlock  of  orthogonality  correction, 

of  FOV  A  the  rate  FOV  A  axes  time  delays, 

in  inertial  measure-  to  the  gyro 

space  ments  frame  other  gyro  error  source 

parameters 


Consistent  with  the  fact  that  the  estimation  algorithm  will  be 
designed  to  prevent  long  term  accumulation  of  rate  errors  and  recognizing 
that  effects  of  deleting  all  gyro  parameters  (except  the  biases)  will 
result  in  most  of  the  errors  being  absorbed  (due  to  high  correlation) 
into  slightly  incorrect  aliased  values  for  the  rate  biases  (which  aren't, 
in  themselves,  of  crucial  interest),  we  have  adopted  the  six  dimensional 


state  vector 


X‘  ■  b.,b,,b,}  (A. 2) 

an  an  an  12  3 

Our  simulations  support  the  conclusion  that  setting  all  other  gyro 
parameters  and  interlocks  to  nominal  apriori  values  does  not  significantly 
degrade  the  Kalman  filter  update*  of  (4'an»9an*4'an)  although  the  forward 
propagations  of  motion  [between  Process  B  outputs]  is  degraded.  For 
the  Phase  I  computer  software  development,  we  have  used  exclusively  the 
six  element  state  vector  of  Eq.  (A. 2). 

The  angular  velocity  measurements  are  considered  to  be  composed  of 


three  terms 


u>1(t) 

ui(t)  + 

bi  +  i  “  1*2.3 

(A. 3) 
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The  raw  signal  w^(t)* **  are  integrated  numerically  into  angular  motion  by 
solving  the  kinematical  differential  equations 
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(4.4)** 


using  4-cycle  Runge-Kutta.  The  equations  (4.4)  are  exact  kinematical 

relationships,  errors  in  [ (t),  6  (t),  4/  ( t )  ]  can  result  only  from  one 

an  an  an 

of  the  following  sources 

(a)  initial  condition  errors 

(b)  errors  in  ($  ,0  , ip  )  estimates 

ag  ag  ag 

(c)  errors  in 

(d)  truncation  and  arithmetic  errors  in  the  numerical  integration 
process. 

Error  sources  (b)  and  (c)  are  correlated  and  correctable  (in  the  Kalman 

filter  algorithm),  errors  (d)  are  not  a  significant  practical  difficulty 

in  the  present  case  except  near  the  geometric  singularity  [ | ©  |  -  90° ] 

an 

of  the  1-2-3  Euler  angle  sequence,  and  since  the  initial  conditions  are 
the  <_  5  arc  second  error  output  of  Process  B,  the  errors  (a)  are  not 
damaging  over  short  time  intervals. 


* 

The  solutions  for  (4>  (t)f0  (t),^  (t)]  are  significantly  improved 

by  judicious  presmoothing  o?  u^(t^?  but  tfte  degree  of  improvement  is  not 
dramatic  over  relatively  short  integration  intervals.  We  have  found  that 
integrating  "raw"  rate  data  (corrected  only  for  biases)  by  four  cycle  Runge- 
Kutta  produces  satisfactory  results  for  the  present  application  (assuming 
about  1  sec. /sec.  noise  for  the  gyro  angular  velocity  measurements). 

** 

We  have  considered  several  other  angular  state  variables  and 
associated  kinematics  (see  section  2.0).  During  Phase  II,  we  will 
conduct  sufficient  studies  on  the  HP  9845S  to  allow  a  final  selection 
to  be  made.  Our  simulations  to  date  have  considered  primarily  the  above 
case. 
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4.2  The  State  Differential  Equations 
STATE  VECTOR; 
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OBSERVATION  VECTOR: 
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lo¬ 


an 


DIFFERENTIAL  EQUATIONS; 
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whare 
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The  two  nontrivial  aubmatrlces  of  partial  derivatives  are  obtained  by 
solving  the  two  differential  equations 

jrl4’ii<t*tk)  1  "  ,F(t)1  ‘♦ii(t*V1 


( 4 . 8d ) 
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and  the  columns  of  (F(t)|  are  given  by  the  following  equations 
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The  Indicated  partlals  In  (4 . 8f )  wore  expllcitely  taken  by  differentiating 
equation  (A. 7a). 


A. 3  The  Kalman  Estimation  Algorithm 

The  classical  extended  Kalman  filter  equations  (see  ref.  A)  in  the 

form 

state JJHDATK  RECURSION 

X^CM-l)  -  X^(k+1)  +  K(k+l)tY(k+l)  -  Y(k+1))  (A. 9) 
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COVARIANCE  AND  KALMAN  GAIN  UPDATE  RECURSIONS 

PR(k+l)  -  »(k+l.k)  Pk(k)»T(k+l,k)  +  Q(k) 
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(4.10) 

(4.11) 


(4.12) 


(4.13) 


Forward  propagated  (apriorl  estimates) 

X^k+l)  ■  X  at  time  tk+^  based  upon  k  data  subsets 

-  forward  Runge-Kutta  integration  equations  in  sec.  4.2 
of  the  previous  estimate  X^k). 

j^k+l* 

Y(k+1)  "  /o (tk+1 )  ■  first  three  elements  of  Xk(tk+.) 


COVARIANCE/PROCESS  NOISE  MATRICES 
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3x3  covariance  matrix  of  errors  associated  with 


($  . 0  .'i'  )»  from  Process  B  output, 

an  an  an 


Q(t)  ■  Process  noise,  to  reflect  uncertainty  in  (4>  ,0  ) 

an  an  an 

arising  due  to  errors  in  the  rate  gyros. 

Our  simulations  indicate  that  angular  rate  errors  (with  Gaussian  white 

2 

noise  and  variance  of  °qyR0  )  result  in  the  following  process  noise  matrix 
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4.4  Some  observations 

The  above  algorithm  is  not  put  forth  as  the  optimal,  final  word  on 
the  subject.  However,  the  simulations  of  section  5  support  the  conclusion 
that  it  is  (in  most  respects)  entirely  satisfactory.  The  primary  shortcoming 


is  the  presence  of  twice  per  orbit  singularities  in  the  1-2-3  angles 
(4  .6  ,<(>)•  As  is  discussed  in  section  2,  this  difficulty  can  be 
readily  overcome  in  a  number  of  ways;  the  best  way  to  overcome  this 
relatively  minox  problem  will  be  established  during  Phase  II. 


S.O  SIMULATIONS 


In  order  to  test  the  logic  and  algorithms  for  all  three  processes, 
it  is  necessary  to  generate  "truth  models"  to  which  simulations  can  be 
compared.  These  truth  models  must  accurately  represent  the  behavior  of 
the  Instrument  with  the  capability  to  simulate  "perfect"  data  as  well  as 
data  with  realistic  (but  known)  error  sources. 

5.1  Process  A  Simulations 

For  Process  A,  we  assume  that  the  star  images  from  the  star  sensors 
results  from  "centroiding"  of  star  energy  images  by  a  70  mm  focal  length 
lens  onto  a  CCD  array,  displaced  slightly  behind  the  image  plane  (2.A25  mm)  to 
allow  for  image  defocus.  Since  the  image  plane  size  is  small  compared 
to  the  focal  length,  we  do  not  include  any  image  plane  distortions  (although 
this  assumption  can  be  readily  modified  for  specific  applications).  Pre¬ 
sumably,  if  there  are  distortions  in  the  final  instrument,  their  effect  will 
be  small  and  assumed  constant  in  time  and,  thus,  will  be  correctable  by 
appropriate  apriorl  calibrations  and  included  in  the  software  for  Process 
A.  We  assume  that  the  true  stellar  image  coordinates  are  given  by  the 
stellar  colinearity  equations  (see  Sec.  3.1). 

For  Process  A  simulations,  we  used  typical  CCD  parameters  for  dark 
current,  sensitivity  and  size.  Assuming  that  a  defocused  star  image  illumi¬ 
nates  a  A  x  A  submatrix  of  pixels,  we  interpolated  the  response  to  determine 
the  image  centroid.  The  number  of  test  cases  was  not  large  but  gave  us 
confidence  that  the  centroid  could  routinely  be  found  with  an  accuracy  to 
about  0.1  pixel.  For  a  pixel  size  of  '0.030  mm,  we  assumed  that  the 
standard  deviation  of  image  errors  is  o'0.003A  mm  (.003A  mm  corresponds 
to  "10  arcseconds  with  a  70  mm  focal  length  lens).  See  Appendix  2  for 
additional  Information  on  parameters  assumed  for  the  CCD.  On  the  basis 
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of  these  tests.  Process  A  output  Is  simulated  by  computing  true  image 
coordinates  and  perturbing  them  slightly  with  simulated  centroiding/ 
measurement  errors  (Gaussian  random  numbers)  of  standard  deviation  .0034 


5.2  Process  B  Simulations 

The  simulations  of  Process  B  consisted  of  determining  the  accuracy  to 
which  Process  B  software  recovered  the  vehicle  orientation  from  a  wide 
variety  of  simulated  star  measurements  from  Process  A.  As  examples,  we 
summarize  ten  cases  in  the  following  discussion.  The  true  orientations 
were  selected  at  random  and  true  stellar  image  coordinates  of  up  to  10 
stars  per  field  of  view  were  perturbed  by  Gaussian  measurement  noise  and 
used  as  input  to  Process  B.  The  initial  angle  estimates  for  vehicle 
orientation  were  displaced  1°  from  the  truth  in  each  angle.  With  this 
initial  estimate,  the  first  cataloged  pair  (after  sorting)  was  usually 
matched  with  the  correct  stored  measured  pair.  Typically,  4-5  least- 
squares  Iterations  were  required  to  bring  the  image  coordinates  of  the 
first  identified  star  pair  into  agreement  with  the  measured  star  pair 
coordinates.  Once  this  was  done,  there  was  usually  no  difficulty  in 
finding  additional  stars  to  confirm  the  match  (i.e.,  the  orientation 
determined  by  a  correctly  identified  star  pair  usually  causes  other 
projected  images  to  nearly  overlay  the  corresponding  measured  images). 

In  several  cases  there  were  fewer  than  expected  confirmations.  This 
occurs  when  one  has  a  particularly  large  error  sample  on  either  the 
first  identified  star  pair  or  on  the  remaining  stars. 

Often,  the  particular  noise  sample  for  image  displacements  can 
impart  a  significant  rotation  of  the  FOV  about  the  boresight.  This  was 
not  a  problem  in  the  combined  solution,  however,  since  the  angle  around 
the  boresight  is  determined  accurately  by  having  the  two  FOV  90°  apart. 
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However,  in  solving  for  the  final  Interlock  angles,  the  rotation  about 
the  boresight  for  FOV(B)  is  often  very  poorly  determined  due  to  this 
noise/geometric  effect.  Since  one  least-  squares  reduction. does  not 
reliably  recover  the  interlock  angles  to  the  desired  precision,  it  is 
necessary  to  average  over  back  least-squares  estimates  of  the  three 
interlocks  in  order  to  smooth  out  these  variations  and  monitor  the  true 
variations  of  interlock  angles.  This  averaging  is  justified  over  up  to 
IS  minutes  of  real  time,  since  any  real  variations  of  interlocks  will  be 
associated  with  an  orbit  period  (thermal  cycling). 

Table  5-1  lists  data  for  the  10  simulation  trials.  For  each  trial, 
the  number  of  stars  used  in  the  least  squares  correction  is  given.  Five 
stars  is  the  current  maximum  number  used  in  the  least-squares  correction 
per  FOV.  Data  for  four  correction  stages  are  give:  1)  solution  of 
FOV(A),  2)  solution  of  FOV(B),  3)  combined  solution  using  FOV(A)  and 
FOV(B),  and  A)  adjustment  of  interlocks  for  FOV(B).  The  deviations 
given  are  deviations  from  the  truth  in  arcseconds.  As  mentioned  earlier, 
the  rotation  about  the  boresight  is  the  most  poorly  determined  angle 
(♦  in  stage  1  and  in  stage  A).  The  interlock  values  given  are 

Sri  Dfl 

single  least-squares  reductions;  these  errors  are  reduced  in  practice, 
proportional  to  the  square  root  of  the  number  of  back  estimates  averaged. 
Note  that  in  nearly  all  combined  solutions  the  deviation  from  the  true 
vehicle  orientation  is  <  10  arcseconds.  Note  the  very  poor  results  for 
example  A.  The  inacuracies  here  are  due  primarily  to  the  fact  that 

91° — near  the  singularity  for  the  1-2-3  Euler  angle  sequence  used  for 
the  Process  B  algorithm.  This  should  not  cause  any  long  term  problems 
since  the  vehicle  will  soon  rotate  to  new  orientations.  Also  in  practice, 
one  will  Introduce  the  "orbiting  frame"  so  that  the  "ao"  angles  will  be 
estimated;  for  this  case,  one  would  generally  be  estimating  small  1-2-3 
angles  (which  will  not  encounter  the  poor  geometry  of  example  A). 


S.3  Process  C  Simulations 

The  primary  function  of  Process  C  that  must  be  tested  is  the  Kalman 
filter  algorithm.  That  is,  we  must  determine  how  well  the  filter  recovers 
the  biases  and  updates  the  optimal  vehicle  orientation.  To  test  this, 
we  computed  and  saved  the  "true"  state,  consisting  of  three  orientation 
angles  for  20  one  minute  intervals.  The  motion  history  is  given  by  the 
rotation  around  the  vehicle  "2"  axis  with  a  rate  of  one  rotation  every 
1.5  hour  plus  a  low  amplitude  sinusoid: 

w  -  w  (1  +  .1  sin((t  -  t  )ut  ))• 
to  o  o 

The  rotation  about  the  other  two  axes  is  nominally  zero.  We  consider  two 
cases: 

(♦,0, -  (0°,  0°,  0°)  and  (♦,0,*)  -  (0°,  91°,  0°). 

an(t-0)  an(t-0) 

For  the  simulations  we  start  with  the  initial  state  estimate  displaced 

1*  from  the  truth  in  all  angles  and  biases  set  to  zero.  Process  B  output  for 
each  interval,  consisting  of  three  discrete  angles,  is  simulated  by 
perturbing  the  true  orientation  angles  by  Gaussian  noise  having  a  standard 
deviation  of  five  arcseconds  [(note,  we  do  not  yet  have  Process  B  and  C 
implemented  in  parallel)].  A  3  x  3  observation  covariance  matrix  associ¬ 
ated  with  the  three  orientation  angles  is  consistently  defined  using  the 
same  standard  deviation.  The  true  angular  rates  are  computed  each  time 
angular  rates  are  needed  by  the  Runge-Kutta  integration  and  the  fixed 
biases  plus  a  random  measurement  noise  was  added  to  all  three  rates. 

For  simplicity,  we  use  identical  "true"  bias  values  for  all  three  rate 
simulations.  The  random  noise  added  to  the  rates  was  uncorrelated  white 
noise  with  a  standard  deviation  of  one  arcsecond  per  second.  The  esti¬ 
mated  biases  from  previous  Kalman  filter  updates  are  subtracted  from  the 
simulated  gyro  rates  described  above  in  each  integration  of  the  kinematic 
equations. 


42 


Figures  5-1  and  5-2  show  the  results  of  two  simulations  typical  of 

those  discussed  above.  Plotted  in  each  figure  is  the  total  deviation 
222  1/2 

& ¥  ■  (A$  +  A0  +  Ai p)  of  the  integrated  state  and  the  filtered  state 

from  the  true  orientation  and  the  estimated  standard  deviation  of  &¥ 

(from  the  updated  covariance  matrix  of  the  Kalman  filter).  The  first  several 
integrated  states  lie  off  the  scale  in  both  figures  because  of  the  inaccurate 
initial  state  estimate  (primarily  the  very  poor  bias  estimates).  Note  that 
after  the  second  and  third  state  updates,  respectively,  that  the  biases  are 
recovered  sufficiently  well  that  each  new  integrated  state  is  displaced  less 
than  30  arcseconds  from  the  true  orientation.  Another  important  point 
to  note  is  that  the  optimal  state  estimate  for  the  orientation  nearly  equals 
the  Process  B  orientation.  This  is  because  the  process  B  angles  have  a  much 
higher  weight,  in  general,  than  the  forward  integrated  angles  from  the 
Runge-Kutta  Integration  of  the  noise-corrupted  gyro  data.  In  fact,  on 
the  scale  of  the  plots  of  Fig.  5-1  and  5-2  the  optimal  state  is  nearly 
Indistinguishable  from  the  Process  B  output  state. 

These  figures  demonstrate  that  the  Kalman  filter  formulation,  in  its 
present  form,  should  provide  sufficiently  good  estimates  for  the  next 
Process  B  cycle.  However,  examination  of  the  recovered  biases  given  by 
the  optimal  estimate  shows  that  the  sequence  of  estimates  have  a  large 
scatter  about  the  true  value,  although  the  average  value  of  each  bias,  over 
time,  is  very  near  to  the  true  value.  This  is  due  primarily  to  the 
similarity  in  magnitude  between  the  baises  and  the  unfiltered  gyro 
measurement  noise.  The  effect  of  the  poorly  determined  biases  is  to  make 
the  next  integrated  state  deviate  more  from  the  truth  than  would  be 
the  case  with  more  precise  bias  values.  As  mentioned  above,  all  steady 
state  deviations  are  less  than  30  arcseconds,  which  is  sufficiently  accurate 
for  re-initiating  Process  B.  For  some  applications,  however,  it  may  be 
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desirable  to  achieve  higher  precision.  In  these  cases,  we  suggest 
either  one,  or  both,  of  two  simple  techniques  to  reduce  errors  due  to 
gyro  errors.  The  first,  and  perhaps  simplest,  is  to  average  the  new 
bias  values  with  several  from  previous  Kalman  filter  solutions  and  use 
the  averages  for  the  next  forward  integration.  Since  the  biases  are 
assumed  constant,  or  at  least  slowly  varying,  this  averaging  should  not 
introduce  significant  error  on  the  time  scale  of  tens  of  minutes.  The 
second  method  is  to  pre-filter  the  noisy  gyro  rate  data  to  smooth  out 
the  worst  of  the  gyro  noise.  Typically,  many  gyro  measurements  should 
become  available  during  each  Process  B  and  C  cycle  (>10  readings  per 
second) .  These  rates  could  be  averaged  mathematically  or  filtered 
electronically  before  being  used  by  the  Runge-Kutta  integration  package. 
Although  these  techniques  will  reduce  the  error,  the  integrated  state 
will  still  be  in  error  due  to  the  errors  in  Process  B  angles  and  thus, 
in  the  optimal  state  which  is  used  as  the  initial  state  for  the  next 
integration.  The  selection  of  any  process  to  reduce  the  integrated  rate 
errors  must  take  into  consideration  not  only  its  effectiveness,  but 
also  its  cost  (mainly,  elapsed  time),  when  introducing  the  algorithm 
into  Process  C.  Significant  time  delays  introduced  may  prove  more 
harmful  (in  a  real- time-mode)  than  the  beneficial  effects  of  better 
angular  rate  knowledge  between  Process  B  outputs.  We  will  address  this 
issue  during  Phase  II. 


6.0  CONCLUDING  REMARKS  AND  OUTLOOK  FOR  PHASES  II  AND  III 
The  following  major  results  have  been  achieved  during  the  Phase  I 
effort: 

A.  A  parallel  process  approach  to  the  real-time  spacecraft 
attitude  estimation  problem  has  been  defined. 

B.  Preliminary  FORTRAN  software  has  been  developed  for: 

1)  Process  B  (star  identification/spacecraft  attitude 

estimation) 

2)  Process  C  (sequential  estimation/motion  integration) 
Simulations  with  these  algorithms  support  the  validity  of 
the  basic  concepts. 

C.  Conversion  of  the  FORTRAN  software  to  run  in  a  microcomputer 
system  (the  HP  9845S)  has  been  initiated,  the  preliminary 
indications  are  most  favorable;  it  appears  that  the  cycle  time 
for  Process  C  (to  process  orientation  angle  estimates  and 
integrate  rate  gyro  data  to  angle  estimates  for  60  real-time- 
seconds)  is  about  10  seconds.  Thus  Process  C  is  presently 
running  about  6  times  as  fast  as  real  time,  using  a  BASIC 
interpreter  -  an  assembly  language  and  optimized  version  of 
this  code  will  run  significantly  faster.  Thus,  there  is  basis 
for  optimism  that  the  real-time  attitude  estimation  objectives 
can  be  achieved,  pending  similar  results  for  the  (much  more 
elaborate)  algorithms  of  Process  B. 

During  Phase  II,  the  following  major  tasks  will  be  addressed: 

A.  Development  of  refined  logic  and  algorithms  for  star  identification 
and  attitude  determination;  the  emphasis  being  to  reduce  the 
probability  of  false  identification  to  a  practically  impossible 
event  (make  the  probability  less  than  10  ^!)  and,  at  the 


same  time,  optimize  the  algorithms  to  facilitate  implementation 
on  the  HP  9845S  (which  will  be  accomplished  during  Phase  II). 

B.  Development  of  improved  algorithms  for  Process  C: 

(1)  Eliminate  geometric  singularities  by  re-formulating  the 
Kalman  filter  in  terms  of  Euler  parameters. 

(2)  Studies  on  various  pre-processors  for  smoothing  the  rate 
gyro  data  under  the  constraints  of  near- real- time,  on¬ 
board  calculations. 

C.  In  preparation  for  interfacing  Processes  B  and  C  in  parallel 
microcomputers  during  Phase  III,  extensive  timing  studies 

and  error  analyses  will  be  conducted  to  assess  the  performance 
of  the  sub-systems  in  the  presence  of  large  a priori  uncertainties 
spurious  images,  and  sharp  attitude  maneuvers.  For  example, 
if  Process  B  can  be  made  to  cycle  several  times  per  minute 
(in  the  steady  state),  then  it  may  be  possible  to  simplify 
the  Process  C  algorithms  without  degrading  the  overall  system 
performance  (since  the  rate  data  need  be  integrated  only  for 
a  few  seconds,  sophisticated  pre-smoothing  may  eliminate 
insignificant  errors  and  thereby  serve  only  to  slow  down 
Process  C’s  cycle  time). 

D.  Study  of  the  evolving  CCD  and  related  sensor  technology  to 
re-evaluate  and  refine  the  sensor  models  and  error  assumptions 
implicit  in  this  research. 

During  Phase  III,  Processes  B  and  C  will  be  implemented  in  parallel, 
intermittently  communicating  microcomputers  and  will  employ  a  "driver" 
simulation  of  Process  A  to  input  simulated  image  centroids,  magnitudes 
and  times.  A  "truth  model"  will  be  developed  and  employed  to  evaluate 
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th«  systems  performance.  The  entire  system  can  than  ba  atudiad  and 
optimised  in  a  raallatic  faahlon.  The  objective  being,  of  course, 
to  gain  the  inaights  necessary  to  design  a  truly  autonomous  and  highly 
reliable  "attitude  navigation"  system.  A  significant  by-product  of 
this  research  program  will  be  a  set  of  well-tested  and  optimised  algorithms 
which  serve  to  automate  ground-based  star  pattern  recognition  and  estimation 
of  spacecraft  orientation  history  (from  star  sensor  and  rate  gyro 
measurements) . 
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APPENDIX  1 


STAR  DATA  BASE  AND  MISSION  CATALOG  CREATION 

The  star  catalog  data  base  system  SKYMAP  (Ref.  5)  has  been  selected 
as  the  master  star  data  base.  The  SKYMAP  catalog  was  developed  from  the 
SAO  catalog  and  other  sources  specifically  for  attitude  determination 
programs  by  NASA-GSFC.  It  is  complete  to  the  eighth  magnitude  in  either 
the  blue  (B)  or  visual  (V)  magnitudes.  Additionally,  the  catalog  contains 
right  ascensions,  declinations,  and,  when  known,  the  spectral  type, 
luminosity  class,  and  amount  of  interstellar  absorption  in  the  V 
wavelength  range. 

The  on-board  (or  mission)  star  catalog  is  divided  into  celestial 
sphere  cells  so  as  to  permit  efficient  microcomputer  access  during  the 
pattern  recognition  process.  In  order  to  keep  storage  requirements  for 
the  mission  catalog  to  a  minimum,  the  cells  do  not  overlap.  The 
placement  of  the  cell  centers  is  given  by  the  polar  angle  0  and  longitude 
X  according  to 

0^  «  cos  *Un)  n  -  0,1,2. ..N  (1.1) 


and 


nj 


ill 

2n+l 


j  "  0,1,2,... 2n 


(-1)"  cos(^),  n  -  0,1,2. ..N 


(1.2) 

(1.3) 


These  formulae  yield  (N+l)  points:  N+l  polar  angles  or  declination 
zones  with  spacing  2«/(2N+l),  and  (2n+l)  equally  spaced  regions  in  each 


zone. 


The  choice  of  N  is  somewhat  arbitrary.  A  large  N  yields  small  cells 
which  would  require  more  than  one  cell  to  be  accessed;  a  small  N  yields 
large  cells  which  would  increase  the  number  of  trials  in  the  pattern 
recognition  process  as  well  as  causing  a  possible  storage  problem.  Taking 


Into  account  the  7®  x  9®  field-of-view,  a  value  of  N  ■  15  was  chosen, 
yielding  256  cells  of  11.6*  x  11.6*  or  larger. 


To  facilitate  computer  access,  the  cells  are  ordered  within  memory 

2 

according  to  a  parameter  n  +j;  a  table  lists  the  starting  relative 
address  of  each  cell  and  the  number  of  stars  in  each.  Thus,  given  a 
boresight  estimate  (6, A),  the  primary  cell  location  is  given  by 
n  -  2[e/ae  +  0.5]  (e<90®) 


(1.4) 

-  2N  +  1  -  2[0/A0  +  0.5]  (6>90°) 

j  -  [A/AA  +  0.5],  (1.5) 


where  [x]  indicates  Integer  arithmetic  (truncation  to  next  smallest 
integer).  The  table  of  cells  is  then  consulted  for  identification  of 
the  appropriate  memory  location.  If  the  estimated  boresight  is  near  a 
cell  boundary  and/or  possesses  a  large  uncertainty,  one  or  more  neighboring 
cells  are  accessed. 

The  CCD  is  assumed  (see  Appendix  2)  to  respond  to  stars  of  I  magnitude 
5  or  lower  -  approximately  5400  stars.  If  these  5400  stars  are  assumed 
to  be  distributed  uniformly  over  the  celestial  sphere,  the  star  density  p 
would  be 

m  5400  stars/sphere _ 

p  41,253  square  degrees /sphere  ^ 

4  0.13  stars/square  degree 

For  a  field-of-view  of  7*  x  9®  *  63  square  degrees,  we  would  expect 

(63  square  degrees)  (0.13  stars/square  degree)  *  8.2  stars* 
in  a  field-of-view  (assuming  uniform  density).  To  obtain  a  measure  of 
the  range  of  the  number  of  stars  actually  detectable  per  field-of-view, 
the  boresight  was  randomly  oriented  over  the  entire  celestial  sphere  100 
times.  For  each  trial,  the  mission  catalog  was  consulted  and  the  number 


Due  to  non-uniform  star  population  of  the  celestial  sphere,  this 
number  decreases  to  about  5  at  the  north  galactic  pole. 


of  stars  in  the  field-of-view  recorded.  The  average  number  of  stars  per 
f ield-of-view  was  six;  in  no  case  were  fewer  than  two  stars  in  the  field- 


APPENDIX  2 


CCD  Instrument  Response  Approximations 
As  mentioned  In  Section  1,  the  outputs  of  Process  A  are  the  Interpolated 
centroids  and  instrument  magnitudes  for  each  valid  star  image.  The  purpose 
of  this  Appendix  is  the  discussion  of  the  approximate  techniques  utilized 
in  the  synthesis  of  these  two  outputs. 

2.1  The  Star  Centroids 


The  centroid  location  (x£,yc)  is  given  by 
E (x . ER . . ) 

v  - 


EER 

ji 


ij 


(2.1) 


and 


z(y^RtJ 
1  3i  J 
EE 

ji  ij 


(2.2) 


where  is  the  A/D  converted  response  level  of  the  pixel  located  at 
^xi,yi^  and  the  summations  are  over  the  square  array  of  pixels  illuminated 
by  the  defocused  star  image  (16  or  25  pixels). 

Typical  cell  size  for  a  CCD  is  approximately  0.030  mm  on  a  side. 

A  CCD  placed  behind  a  70  ram  focal  length  lens  (proposed  for  one  CCD  star¬ 
sensor)  gives  a  resolution  of  approximately  1.5  arc-minutes  for  a  focused 
image.  Uhen  spread  over  a  4  x  4  or  5  x  5  cell  pattern,  the  resolution 
with  which  the  centroid  can  be  located  has  been  found  to  be  about  6  arc- 
sec.  For  double  stars.  Process  A  would  produce  image  coordinates  for  a 
single  star  but  with  poorly  determined  image  coordinates  (a  weighted 
mean  of  the  two  stars). 

Detections  of  double  stars  should  not  be  used  in  Process  B  since 
they  would  result  in  poor  orientations.  One  solution  to  this  problem  is 
to  delete  from  the  mission  catalog  all  star  pairs  with  separations  less 
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Chan  some  tolerance  ('6  arc-oln.  In  this  case).  There  are  sufficient 
stars  in  the  catalog  that  this  deletion  should  not  seriously  degrade 
performance.  Some  additional  time  in  Process  B  will  be  used  trying 
unsuccessfully  to  pair  measured  double  stars  with  catalog  stars,  but 
since  detection  of  double  stars  will  be  a  relatively  rare  event,  this 
time  penalty  should  not  be  a  significant  practical  problem. 

2.2  CCD  Magnitude  Response 
2.2.1  Magnitude  Conversion 

Due  to  both  the  different  spectral  qualities  of  various  stars  and  the 
peculiarities  in  the  unfiltered  CCD  response  (the  primary  sensitivity  is 
to  red  or  near  intra-red  radiation)  two  stars  of  the  same  visual  (V) 
magnitude  (for  example)  may  cause  different  CCD  response.  Hence  none 
of  the  cataloged  star  magnitudes  may  be  used  directly.  Rather,  the 
magnitudes  of  the  stars  must  be  properly  transformed  (using  the  spectral 
properties  contained  in  the  master  star  catalog  SKYMAP)  prior  to 
insertion  in  the  mission  catalog.  Simply  stated,  the  mission  catalog 
must  contain  an  "instrument  magnitude"  for  each  star. 

It  has  been  decided  to  convert  the  V  magnitudes  of  SKYMAP  to  I 
magnitudes  and  utilize  an  I  filter  placed  over  the  CCD  array.  The 
following  points  support  this  decision: 

(1)  Given  the  information  in  SKYMAP,  one  could,  in  principle, 
perform  a  magnitude  conversion  from  V  to  a  CCD  magnitude. 
However,  this  would  require  the  choice  of  a  particular  CCD 
detector  in  order  to  determine  its  response  characteristics  in 
the  laboratory. 


(2)  The  response  functions  of  typical  CCD's  are  quite  broad,  a 
fact  which  makes  a  rigorous  conversion  to  a  CCD  magnitude 


difficult'  la  light  of  the  complicated  stellar  spectral  features 
In  the  blue  wavelength  region.  A  detailed  description  of  the 
spectra  would  be  requled.  The  I  filter,  on  the  other  hand, 
is  confined  to  the  red  wavelengths  where  the  star  spectra 
are  relatively  snooth. 

(3)  The  1  filter  response  peak  is  near  that  of  typical  CCD’s. 

In  addition,  it  overlaps  the  main  peaks  of  the  commonly  accepted 
"typical  CCD  response"  (Ref.  6).  Hence,  an  1  filter  placed 
over  the  CCD  array  would  serve  to  limit  the  wings  of  the  CCD 
response  and  still  provide  adequate  through-put  for  sensitivity. 

(A)  Information  exists  for  converting  V  magnitudes  to  I  magnitudes. 
The  transformation  requires  only  spectral  type  and  luminosity 
class  -  both  readily  available  from  SKYMAP. 

In  the  ideal  case,  the  set  of  detectable  stars  exactly  matches  the 
catalog.  Since  this  is  not  possible,  it  Is  desirable  to  maximize  the 
completeness  of  the  catalog  to  some  rather  faint  magnitude  to  insure  that 
most  detectable  stars  are  contained.  It  is  important  to  note  that  stars 
which  are  faint  in  V  may  be  relatively  brighter  at  red  wavelengths.  As 
will  be  demonstrated  in  the  next  sections,  a  limit  of  magnitude  5  in  1 
seems  to  be  a  reasonable  limit  for  the  CCD  configuration  assumed  for  the 
present  study.  The  8th  magnitude  limit  of  SKYMAP  assures  that  the 
mission  catalog  listing  will  be  sufficient. 


2.2.2  The  I  Magnitude  Conversion  Method 


Two  external  items  of  information  are  needed  for  the  magnitude 
conversion.  Values  of  V^g  "  Igbg  were  obtained  from  Johnson  [7],  whose 
table  contains  listings  for  three  luminosity  classes  (I,  III,  &  V  -  super 
giants,  giants  and  main  sequence)  and  extends  over  most  spectral  types. 
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The  second  value  needed  is  a^/a^,  the  ratios  of  absorption  in  I  to 
absorption  in  V,  expressed  in  magnitude  and  as  a  function  of  spectral 
type  by 

ax  log10[(/o”l(A)E(A)(l  -  I<A)dA)/(/o“l(A)E(A)dA>]  (2 

®V  log10[(/o*V(A)E(A)(l  -  I(A)dA)/(/o"v(A)E(A)dA)] 

Where  1(A)  and  V(A)  are  the  filter  responses,  E(A)  is  the  star  energy 

function  and  1(A)  is  the  relative  absorption  function  and  A  is  the  wavelength. 

(see  SKYMAP  discription  [5]  for  details).  As  pointed  out  in  SKYMAP,  this 

ratio  is  nearly  constant  over  spectral  type  for  narrow  or  intermediate  band 

filters.  To  calculate  this  ratio  the  Planck  energy  function  was  used  to 

model  the  stellar  flux.  Although  this  is  not  precisely  valid,  forming 

the  ratio  should  lead  to  quite  accurate  results.  The  temperatures  used 

were  those  given  by  Johnson  [7]  and  no  distinction  was  made  by  luminosity 

class.  The  values  for  absorption  were  taken  from  Figure  3.2  in  the 

SKYMAP  description  by  assuming  absorption  in  magnitude  is  a  linear 

•  • 

function  of  wavelength  over  the  range  of  interest  (4800  A  -  10000  A) , 

a(A)  -  1.77  x  10"A  x  A (A)  +  1.77,  (2.4) 

Where  a(A)  is  absorption  in  magnitude  at  wavelength  A.  The  results  of 
these  calculations  were  that  a^/a^  varied  from  0.25  to  0.32.  The  same 
value  of  was  used  for  I,  III  and  V  luminosity  class  stars  at  a 

given  spectral  type. 

2.2.3  SKYMAP  DATA 

The  data  from  SKYMAP  needed  for  the  magnitude  conversion  consists 
of  apparent  visual  magnitude,  spectral  type,  luminosity  class  and 
absorption  in  V.  These  data,  with  the  exception  of  av,  are  given  for 


difficult  in  light  of  the  complicated  stellar  spectral  features 
in  the  blue  wavelength  region.  A  detailed  description  of  the 


spectra  would  be  requled.  The  I  filter,  on  the  other  hand, 
is  confined  to  the  red  wavelengths  where  the  star  spectra 
are  relatively  smooth. 

(3)  The  I  filter  response  peak  is  near  that  of  typical  CCD's. 

In  addition,  it  overlaps  the  main  peaks  of  the  commonly  accepted 
"typical  CCD  response"  (Ref.  6).  Hence,  an  I  filter  placed 
over  the  CCD  array  would  serve  to  limit  the  wings  of  the  CCD 
response  and  still  provide  adequate  through-put  for  sensitivity. 

(4)  Information  exists  for  converting  V  magnitudes  to  I  magnitudes. 
The  transformation  requires  only  spectral  type  and  luminosity 
class  -  both  readily  available  from  SKYMAP. 

In  the  ideal  case,  the  set  of  detectable  stars  exactly  matches  the 

catalog.  Since  this  is  not  possible,  it  is  desirable  to  maximize  the 

completeness  of  the  catalog  to  some  rather  faint  magnitude  to  Insure  that 

most  detectable  stars  are  contained.  It  is  Important  to  note  that  stars 

which  are  faint  in  V  may  be  relatively  brighter  at  red  wavelengths.  As 

will  be  demonstrated  in  the  next  sections,  a  limit  of  magnitude  5  in  I 

seems  to  be  a  reasonable  limit  for  the  CCD  configuration  assumed  for  the 

present  study.  The  8th  magnitude  limit  of  SKYMAP  assures  that  the 

mission  catalog  listing  will  be  sufficient. 

2.2.2  The  I  Magnitude  Conversion  Method 

TWo  external  items  of  information  are  needed  for  the  magnitude 

conversion.  Values  of  V  .  -  I  .  were  obtained  from  Johnson  [7],  whose 

abs  abs 

table  contains  listings  for  three  luminosity  classes  (I,  III,  &  V  -  super 
giants,  giants  and  main  sequence)  and  extends  over  most  spectral  types. 
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most  scars.  It  was  necessary  Co  collapse  catagorses  of  luminosity  class 
and  some  spectral  types  since  the  table  from  Johnson  is  limited.  Collapsing 
is  justified  in  most  cases  because  either  the  category  contains  few  stars 
and/or  the  properties  are  similar  to  those  of  listed  star  types.  The 
following  combinations  were  used: 


Luminosity  class  III,  spectral  type  M 
(all  are  variants  of  M  III) 


Lum  class  III,  spectral  type  09 

(all  are  hot  giant  stars  similar  to  09) 


0 
H 


^a’ab’b  ) 

No  Luminosity  class  ■*  V  (most  stars  are  V  stars) 

No  subinterval  in  spectral  type  -*-5  (i.e.,  A  becomes  A3) 

No  absorption  given  set  to  0 
No  spectral  type  exclude 
No  my  -*•  exclude 

Given  the  spectral  type  and  luminosity  class,  the  value  of  -  l>bg  and  a^/ay 


were  found  by  interpolation  in  the  table.  Then: 

*1  “  "V  “  *Vaba  “  ^bs*  ‘  *V^  -  *I,*V** 
where  my  and  ay  are  from  SKYMAP. 


(2.5) 
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In  order  to  establish  a  reasonable  magnitude  limit  for  a  CCD  sensor 
we  1)  determine  the  flux  in  the  filter  bandpass  for  some  standard  star 
at  the  earth's  atmosphere,  and  2)  multiply  by  appropriate  factors  dictated 
by  the  sensor. 

The  most  direct  way  to  obtain  a  flux  estimate  would  be  to  observe 

known  stars  from  space  with  the  CCD  sensor.  Barring  this,  ground-based 

observations  of  stars  with  varying  zenith  angles  could  yield  flux 

estimates  outside  the  earth's  atmosphere. 

Our  approximate  method  was  to  numerically  integrate  the  surface 

flux  distribution  of  a  K7  V  star  model  atmosphere  over  the  I  bandpass. 

The  absolute  I  magnitude  of  such  a  star  is  approximately  6.2  (V  »  8.1, 

V  -  I  «  1.93).  The  radius  is  given  by  log  R*/R  ■  -0.11  where  R  » 

sun  R  sun 

6.96  x  lO^cm.  The  surface  flux  is  scaled  by  -  x  ^  «  3.08  x  10 

K  JLU 

sun  pc 

where  1  pc  “  3*08  x  10*  cm.  The  result  of  integrating  and  scaling  is: 

2 

'1024  photons /cm  sec. 

Typical  scale  factors  are: 

2 

-Lens  area:  26.7  cm 

-CCD  response  peak  efficiency  ■  ,60 

-1  filter  transmission  peak  efficiency  »  0.85 

-CCD  effective  area  utilization  “  0.46 
(based  on  Fairchild  488  x  380  cell  CCD) 

-Integration  time  *  0.1  sec. 

If  we  desire  a  minimum  of  7500  photons/star  for  a  sufficient  signal-to- 
noise  ratio,  we  compute  the  1  magnitude  limit  of: 

“limit  “  6‘2  ♦  2’5  lo*  755o 


Mote  Chat  many  factors  are  uncertain  or  could  be  altered.  Integration 
time  could  be  increased  to  1  sec  to  give  a  limit  of  6.0.  We  have  chosen 
5.0  to  be  the  cutoff  magnitude  since  this  seems  obtainable  fcnd  gives 
approximately  5400  stars,  a  sufficient  number  for  the  pattern  recognition 
process  to  work  reliably. 

We  also  note  that  model  atmospheres  for  a  variety  of  spectral  types 
could  be  used  to  repeat  the  above  calculation  to  yield  a  more  precise 
magnitude  limit. 

The  magnitude  limit  is  flexible  since  the  integration  time  for  the 
star  sensor  is  variable  over  a  wide  range.  If  the  integration  time  is 
changed  by  a  factor  of  10  the  magnitude  limit  is  changed  by  2.5  mag. 

In  addition,  the  dynamic  range  of  typical  CCD  arrays  is  200  or  about  6 
magnitudes.  The  response  is  linear  over  the  range  to  allow  accurate 
magnitude  calibration  and  detection. 
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